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Abstract: Aeolian sandy soil in mining areas exhibits intense evaporation and poor water retention 
capacity. This study was designed to find a suitable biochar application method to improve soil water 
infiltration and minimize soil water evaporation for aeolian sand soil. Using the indoor soil column 
method, we studied the effects of three application patterns (A (0—20 cm was a mixed sample of 
mixed-based biochar and soil), B (0—10 cm was a mixed sample of mixed-based biochar and soil and 
10-20 cm was soil), and C (0—10 cm was soil and 10—20 cm was a mixed sample of mixed-based biochar 
and soil)), four application amounts (0% (control, CK), 1%, 2%, and 4% of mixed-based biochar in dry 
soil), and two particle sizes (0.05—-0.25 mm (S1) and <0.05 mm (S2)) of mixed-based biochar on water 
infiltration and evaporation of aeolian sandy soil. We separately used five infiltration models (the Philip, 
Kostiakov, Horton, USDA-NRCS (United States Department of Agriculture-Natural Resources 
Conservation Service), and Kostiakov-Lewis models) to fit cumulative infiltration and time. Compared 
with CK, the application of mixed-based biochar significantly reduced cumulative soil water infiltration. 
Under application patterns A, B, and C, the higher the application amount and the finer the particle size 
were, the lower the migration speed of the wetting front. With the same application amount, cumulative 
soil water infiltration under application pattern A was the lowest. Taking infiltration for 10 min as an 
example, the reductions of cumulative soil water infiltration under the treatments of A2%g2, A4%1), 
A4% 2), AL% G1, C2%s1, and B1% 61) were higher than 30%, which met the requirements of loess soil 
hydraulic parameters suitable for plant growth. The five infiltration models well fitted the effects of the 
treatments of application pattern C and S1 particle size (R?>0.980), but the R? values of the Horton 
model exceeded 0.990 for all treatments (except for the treatment B2%g2)). Compared with CK, all other 
treatments reduced cumulative soil water infiltration, except for B4%z). With the same application 
amount, cumulative soil water evaporation difference between application patterns A and B was small. 
Treatments of application pattern C and S1 particle size caused a larger reduction in cumulative soil water 
evaporation. The reductions in cumulative soil water evaporation under the treatments of C4% 1, 
C4% 2), C2% 1, and C2%s2) were over 15.00%. Therefore, applying 2% of mixed-based biochar with S1 
particle size to the underlying layer (10—20 cm) could improve soil water infiltration while minimizing soil 
water evaporation. Moreover, application pattern was the main factor affecting soil water infiltration and 
evaporation. Further, there were interactions among the three influencing factors in the infiltration 
process (application amountXparticle size with the most important interaction), while there were no 
interactions among them in the evaporation process. The results of this study could contribute to the 
rational application of mixed-based biochar in aeolian sandy soil and the resource utilization of urban and 
agricultural wastes in mining areas. 


“Corresponding author: JIA Jianli (E-mail: jj1@cumtb.edu.cn) 

Received 2021-08-03; revised 2021-12-29; accepted 2022-01-06 

© Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Science Press and Springer-Verlag GmbH Germany, 
part of Springer Nature 2022 


http://jal.xjegi.com; www.springer.com/40333 


ZOU Yiping et al.: Effects of mixed-based biochar on water infiltration and evaporation... 


Keywords: biochar; water infiltration; water evaporation; aeolian sand soil; mining areas 


Citation: ZOU Yiping, ZHANG Shuyue, SHI Ziyue, ZHOU Huixin, ZHENG Haowei, HU Jiahui, MEI Jing, BAI Lu, JIA 
Jianli. 2022. Effects of mixed-based biochar on water infiltration and evaporation in aeolian sand soil. Journal of Arid 
Land, 14(4): 374-389. https://doi.org/10.1007/s40333-022-0060-6 


1 Introduction 


In most arid and semi-arid regions in the world, soils are characterized by high sand and low 
organic matter contents, which harm the physical and chemical properties of soils. Specifically, 
these soils feature high water infiltration rate, low water holding capacity, large evaporation, and 
low fertility and organic matter content, leading to low water supply and utilization efficiency 
(Alaa et al., 2017). Generally, low soil organic matter content results in poor water holding 
capacity and weak soil structure (Breus et al., 2014). Studies have shown that the physical and 
chemical properties of soils can be improved by adding organic and inorganic amendments, such 
as biochar (Tang et al., 2015; Basanta et al., 2017; Zhou et al., 2020). 

Recently, it has been reported that biochar from pyrolysis of organic waste can be considered 
an alternative amendment that has the characteristics of porous structure, large specific surface 
area, and rich functional groups (Wang et al., 2021a, b), which affect the bulk density, porosity, 
aggregation, and hydraulic properties of soils (Lehmann, 2007; Lehmann et al., 2011; Ahmad et 
al., 2012; Kinney et al., 2012; Li, 2019). Studies have shown that biochar can reduce soil bulk 
density (Laird et al., 2010; Githinji, 2013), increase soil porosity (Tammeorg et al., 2014), change 
the number of soil aggregates and their stabilities (Mukherjee et al., 2014; Zhang et al., 2015), 
and affect soil hydraulic properties (Wang et al., 2015; Bao, 2020; Razzaghi et al., 2020). These 
effects of biochar application, however, depend on various factors, such as application pattern, 
application amount, particle size, type of biochar, and soil type (Qi et al., 2014; Xie et al., 2016; 
Wang et al., 2018; Wang et al., 2019). At present, researchers are focusing on exploring the effects 
of application patterns and application amounts of biochar on soil hydraulic properties; it has been 
found that an appropriate amount of biochar application can improve the infiltration capacities of 
soils, and different application patterns of added biochar have significant impacts on soil water 
infiltration and evaporation (Burrell et al., 2016; Li et al., 2016, 2018; Lim et al., 2016; Liu et al., 
2016). However, the effects of particle sizes of biochar on soil hydraulic characteristics are 
uncertain (Xie et al., 2016; Wang et al., 2019), and most studies focus on particle size above 0.25 
mm, and less studies below 0.25 mm (Wang et al., 2019; Kim et al., 2021). Therefore, the effects 
of application patterns, application amounts, and particle sizes of biochar on soil water infiltration 
and evaporation are not clear. Moreover, in arid and semi-arid regions, it is unclear whether 
suitable application patterns, application amounts, and particle sizes of biochar can significantly 
improve soil water infiltration and effectively reduce water evaporation from aeolian sandy soil, 
especially in mining areas with degraded ecological environments. Exploring this issue is very 
important for determining whether biochar can improve the soil quality of mining areas with 
insufficient water retention capacity. Additionally, the raw materials of biochar are mostly wood 
and straw; the sludge from urban sewage treatment plants usually contains large amounts of 
organic matter and is a suitable raw material for preparing biochar. Thermal conversion of sludge 
through pyrolysis reduces the volume of sludge, kills pathogens and parasites, and realizes the 
safe disposal and resource utilization of sludge (Sun et al., 2018). 

In this study, we used co-pyrolysis of urban domestic sludge and agricultural straw to prepare 
mixed-based biochar. Through a one-dimensional vertical infiltration experiment, we explored the 
effects of mixed-based biochar on the soil water infiltration and evaporation characteristics of 
aeolian sandy soil in a mining area of China. Specifically, we set up three application patterns for 
mixed-based biochar, four application amount gradients, and two particle size ranges. We adopted 
a full-factorial design experiment with a total of 19 treatments. The following goals are expected 
to be achieved: (1) exploring the effects of application patterns, application amounts, and particle 
sizes of mixed-based biochar on soil water infiltration and evaporation; and (2) determining the 
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optimal combination of application pattern, application amount, and particle size of mixed-based 
biochar to improve water infiltration and minimize water evaporation in aeolian sandy soil. The 
results of this study could contribute to the rational application of biochar in aeolian sandy soil 
and the resource utilization of urban and agricultural wastes in mining areas. 


2 Materials and methods 


2.1 Study area 


In this study, we collected aeolian sandy soil from the collapse area of the Daliuta Coal Mine 
(110°26’N, 39°29’E) located in Daliuta Town, Shenmu County, Shaanxi Province, China. The 
mining area belongs to the Shendong mining area, which is the transition zone between the 
southern edge of the Mu Us Sandy Land and the northern edge of the Loess Plateau, with a poor 
environment and fragile ecology. The annual precipitation is between 251.3 and 646.5 mm, and 
the evaporation is as high as 1788.4 mm. The precipitation is highly variable and concentrated in 
summer. The average annual temperature is about 7.0°C. 


2.2 Sample collection 


Soil samples were collected from 14 sampling points in Shendong Central Coal Mine, such as 
Shangwan, Daliuta, and Halagou. We collected aeolian sandy soil samples at depths of 0-50 cm 
in July 2020 by using the biological community sampling method. After removing the impurities, 
the soil samples were naturally dried in the dark and then ground through a 2-mm sieve for 
analysis. A laser particle size analyzer (Malvern 2000, Mettler Toledo, Malvern, USA) was used 
to determine the mechanical components of aeolian sandy soil. The mass fractions of sand, silt, 
and clay were 85.52%, 14.44%, and 0.04%, respectively. Other physical and chemical properties 
of the tested soil are shown in Table 1. 

We took sludge from the Gaobeidian Sewage Treatment Plant in Beijing City, China. After the 
sludge was air-dried at room temperature, it was crushed with a grinder and passed through a 
60-mesh screen. The dried powder after sieving was stored in a dryer to prevent moisture 
adsorption. We collected corn stalks from the suburbs of Beijing City. After 10 d in the sun, the 
corn stalks were chopped with a chopper, passed through a 60-mesh screen, and stored in a 
desiccator. Mixed-based biochar was prepared by mixing sludge and corn stalks with a ratio of 
7:3 (w/w) and carbonizing at 550°C with a residence time of 120 min (Table 2). 


Table 1 Physical and chemical properties of tested aeolian sandy soil 


H EC CEC TC SOC TN AN TP AP Number of 
P (uS/cm) (cmol/kg) (g/kg) (g/kg) (g/kg) (mg/kg) (g/kg) (mg/kg) microorganisms (x10*/g) 
7.77 309.91 5.00 18.70 10.67 0.65 52.68 0.51 32.84 151.99 


Note: EC, electrical conductivity; CEC, cation exchange capacity; TC, total carbon; SOC, soil organic carbon; TN, total nitrogen; AN, 
available nitrogen; TP, total phosphorus; AP, available phosphorus. 


Table 2 Physical and chemical properties of mixed-based biochar 


H CEC BET C H N O P K Zn Cu Cr Average pore Pore volume 
ph (cmol/kg) (m/z) (%) (%) (%) (%) (g/kg) (g/kg) (mg/kg) (mg/kg) (mg/kg) size (nm) (cm*/g) 
8.02 19.12 72.03 61.02 3.51 066 893 4.81 9.50 183.99 8.31 1.79 9.77 0.18 


Note: BET, Brunner-Emmet-Teller; C, carbon; H, hydrogen; N, nitrogen; O, oxygen; P, phosphorus; K, potassium; Zn, zinc; Cu, 
copper; Cr, chromium. 


2.3 Experimental design 


In this study, we set up three application patterns (Fig. 1), four application amount gradients, and 
two particle sizes of mixed-based biochar. We adopted a full-factorial design experiment with a 
total of 19 treatments (Table 3), and each treatment was repeated three times. Application patterns 
of mixed-based biochar were as follows: (1) 0-20 cm was a mixed sample of mixed-based 
biochar and soil (A); (2) 0-10 cm was a mixed sample of mixed-based biochar and soil and 10-20 
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cm was soil (B); and (3) 0-10 cm was soil and 10-20 cm was a mixed sample of mixed-based 
biochar and soil (C). The four application amounts of mixed-based biochar in dry soil were 0% 
(CK, control), 1%, 2%, and 4%. The two particle sizes of mixed-based biochar were 0.05—0.25 
mm (S1) and <0.05 mm (S2). According to the prefilling test, the filling bulk densities of soil and 
mixed-based biochar were 1.56 and 0.40 g/cm, respectively. Since mixed-based biochar reduces 
soil bulk density, we kept the volume of filling soil constant. We determined the method used for 
filling the mixed soil bulk density according to the application amounts of mixed-based biochar, 
soil bulk density, and soil addition ratio. Therefore, soil bulk densities of soil with 1%, 2%, and 
4% mixed-based biochar were 1.52, 1.48, and 1.40 g/cm’, respectively. 


E Mixture of mixed-based biochar and soil E Soil 


Surface layer 
(0-10 cm) 
Underlying layer 
(10-20 cm) 
CK A B C 


Treatment 


Fig. 1 Schematic diagram showing the application patterns of mixed-based biochar. CK, control. Application 
patterns for mixed-based biochar were as follows: (1) 0-20 cm was a mixed sample of mixed-based biochar and 
soil (A); (2) 0-10 cm was a mixed sample of mixed-based biochar and soil and 10-20 cm was soil (B); and (3) 
0-10 cm was soil and 10-20 cm was a mixed sample of mixed-based biochar and soil (C). 


Table 3 Full-factorial design experiment used in this study 


Treatment Application pattern Application amount Particle size 
CK 
Al% 1) = S1 
A1%s2) S2 
A2% 51) S1 
A 2% 
A2% 52) S2 
A4% 51) op S1 
A4% 52) S2 
B1% 1) ia S1 
B1%isz S2 
B2% 1) S1 
B 2% 
B2% 2) $2 
B4% 51) iu S1 
B4% 52) $2 
C1%s1) i S1 
C1 % 52) $2 
C2%s1) č ia S1 
C2%s2) $2 
C4% 51) 4a S1 
C4% 52) $2 


Note: CK, control. Application patterns of mixed-based biochar were as follows: (1) 0-20 cm was the mixed sample of 
mixed-based biochar and soil (A); (2) 0-10 cm was the mixed sample of mixed-based biochar and soil and 10-20 cm was soil 
(B); and (3) 0-10 cm was soil and 10—20 cm was the mixed sample of mixed-based biochar and soil (C). $1, 0.05—0.25 mm; S2, 
<0.05 mm. 
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2.4 Determination of soil water infiltration and evaporation 


2.4.1 Preparation of the soil column 


The soil column used in the experiments was a transparent acrylic column with an inner diameter 
of 7.0 cm and a height of 30.0 cm. Three layers of gauze sand and one layer of gravel were placed 
in the bottom of the soil column to prevent the loss of soil particles, and then the soil was added. 
Vaseline was applied evenly and thinly on the pipe wall to reduce pipe wall effects on soil water 
infiltration. Each soil column was compacted every 5.0 cm, and mixed-based biochar and soil 
samples were calculated and weighed separately, mixed evenly, loaded into the soil column, and 
carefully "flashed" at the layered interface to ensure close contact between the soil layers to avoid 
delamination. After installing the soil column, a scraper was used to smooth the surface. 


2.4.2 Determination of soil water infiltration and evaporation 


The one-dimensional constant head vertical ponding infiltration method was used to measure the 
soil water infiltration parameters. The infiltration device was composed of a Mahalanobis flask 
and a soil column (Fig. 2). Water was added to the Mahalanobis flask. A rubber hose was used to 
connect the water inlet of the soil column to the outlet hole of the Mahalanobis flask. We adjusted 
the relative height of the soil column and Mahalanobis flask to make the outlet hole of the 
Mahalanobis flask discharge water. At this time, the head height was 2.5 cm. After the beginning 
of infiltration, we continuously recorded the descending depth of the vertical wetting front and the 
descending height of the water level in the Mahalanobis flask for a certain period. When the 
wetting front moved to the bottom of the soil column, recording was stopped, but the water 
supply continued until the soil column was saturated. Then, all treated soil columns were placed 


Mariotte bottle 


-—— Soil column 


Rubber hose | 


sii ty 
| aeee Water head 


Spout plug —§ o — 


Thickness of soil 


Gravel U- 
Spreader plate z : 
i 


Fig. 2 Demonstration of the device used in determinations of soil water infiltration and evaporation 
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indoors, and an evaporation test was carried out in a relatively stable environment. During this 
period, the average temperature in the control room was 23.4°C, and the average relative 
humidity was 22%. Soil water evaporation losses were measured with an electronic scale at 17:00 
(LST) every afternoon, and soil water evaporation was recorded continuously for 53 d. The 
formula used for calculating daily soil water evaporation is as follows: 


E=M,~x10/(ar’), (1) 


where E is the daily evaporation of the soil column (mm); Mz is the daily mass change of the soil 
column (g); and r is the inner radius of the soil column (cm). 


2.5 Data processing and analysis 


All data for the experiments were the averages of repeated determinations. Excel 2016 was used 
for data processing, Origin 2018 software was used for figure mapping, SPSS 20.0 software was 
applied for the simulations of soil water infiltration parameters and statistical analyses, and least 
significant difference (LSD) method was utilized for significance testing (P<0.05). 


3 Results 


3.1 Effects of application patterns, application amounts, and particle sizes of mixed-based 
biochar on the wetting front 


One-dimensional vertical movement of the wetting fronts in all treatments increased gradually 
with infiltration time, but the times taken to reach the bottom of the soil column were different, 
and the migration distance of the wetting front for each treatment also differed with time. 
Compared with CK, all other treatments (except for C1%s1), which promoted the downward 
migration of the wetting front) slowed down the migration speed of the wetting front. Under 
application patterns A and B of mixed-based biochar, the higher the application amount and the 
finer the particle size were, the lower the migration speed of the wetting front. The degree of the 
migration speed slowed down by mixed-based biochar decreased in the order of 
4% (52)>2%(52)>4%51)> 1 %is2) and 2% s1)>1% s1)>CK. Figure 3 shows that the treatment of 4%,s2) 
had a significantly higher mitigation effect on the downward shift of the wetting front than other 
arbitrary treatments. Under application pattern C, the downward migration speed of the wetting 
front was consistent with those under application patterns A and B, and the descending order of 
the migration speed was 4% s2)>2% s2), 4% s1), and 2% s1)>1%s2) and CK>1% 31). The results 
showed that the treatment of 4%,s2) was still the best combination for slowing down the migration 
speed of the wetting front. In short, under application patterns A, B, and C, the higher the amount 
of mixed-based biochar and the finer the particle size of mixed-based biochar were, the more 
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Infiltration time (min) Infiltration time (min) Infiltration time (min) 
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Fig. 3 Dynamic changes in wetting fronts under different application patterns, amounts, and particle sizes of 
biochar. Application patterns of mixed-based biochar were as follows: (1) 0-20 cm was a mixed sample of 
mixed-based biochar and soil (A); (2) 0-10 cm was a mixed sample of mixed-based biochar and soil and 10—20 
cm was soil (B); and (3) 0-10 cm was soil and 10-20 cm was a mixed sample of mixed-based biochar and soil 
(C). S1, 0.05-0.25 mm; S2, <0.05 mm. The four application amounts of mixed-based biochar in dry soil were 0% 
(control, CK), 1%, 2%, and 4%. 
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pronounced the slowing effect of mixed-based biochar on the migration of the wetting front. With 
the 1% application amount and S1 particle size, application patterns A, B, and C all slowed down 
the migration speed of the wetting front, and the effects decreased in the order A>B>C. However, 
with the 1% application amount and S2 particle size, both application patterns A and B slowed 
down the downward migration of the wetting front and the effects were similar, but application 
pattern C promoted the downward migration of the wetting front. Application patterns A, B, and 
C all slowed down the migration speeds of the wetting fronts with application amounts of 2% and 
4%. In conclusion, for mixed-based biochar with the three application amounts and two particle 
sizes, both the effects of application patterns A and B on the migration of the wetting fronts were 
significantly greater than that of application pattern C. 

Considering application patterns, application amounts, and particle sizes of mixed-based 
biochar, for the infiltration time of 10 min, only the treatment of C1%,s1) promoted soil water 
infiltration. The migration distances of the wetting front under other treatments were all smaller 
than those under CK and reached a significant level (P<0.05), which could increase the soil water 
absorption and slow down the downward migration speed of water in the soil to a certain extent, 
especially for the treatment of A4% s2). 

The relationship between the wetting front and time followed the power function: F=pt" (where 
F is the wetting front, p and v are the empirical parameters, and t is the time). The fitting results in 
Table 4 show that the power function simulated the migration law of the mixed soil moisture peak 
under different application patterns, application amounts, and particle sizes well. For each 
treatment, the determination coefficient (R?) was higher than 0.911, and the test was significant at 
the P<0.05 level. There was no clear rule for each treatment under application patterns A and B, 
but migration was significantly reduced compared with CK, and the treatment of 4% s2) had the 
significant reduction effect. Under application pattern C, the value of p first increased slowly and 
then sharply increased with increases in application amounts. With increases in application 
amounts of mixed-based biochar, the power index v under each treatment (with exception of 
B4%,s2)) showed significant decreasing trends for applications A, B, and C; the smaller the 
particle size was, the lower the value of the power index v. This indicated that application 
patterns, application amounts, and particle sizes of mixed-based biochar significantly impacted 
the initial water infiltration process that dominated by matrix potential, and they all played a vital 
role in attenuating the latter wetting front. 


Table 4 Fitting results for the wetting front and infiltration time under different application patterns, application 
amounts, and particle sizes of mixed-based biochar 


Application pattern A Application pattern B Application pattern C 

Treatment > - m 

p v R p v R p v Rr 
CK 3.867 0.685 0.998 3.867 0.685 0.998 3.867 0.685 0.998 
1% s1) 2.240 0.746 0.994 2.582 0.759 0.991 2.013 1.079 0.951 
1% (52) 2.949 0.602 0.997 2.200 0.683 0.985 2.542 0.928 0.967 
2% s1) 2.108 0.723 0.991 2.878 0.608 0.974 2.023 0.930 0.951 
2% 52) 2.048 0.616 0.999 2.349 0.563 0.988 2.649 0.782 0.938 
4% 51) 2.406 0.594 0.990 2.604 0.565 0.985 3.412 0.685 0.947 
4% (52) 1.851 0.539 0.952 1.442 0.614 0.943 3.553 0.599 0.911 


Note: p and v are the empirical parameters, and R? is the determination coefficient. 


3.2 Effects of application patterns, application amounts, and particle sizes of mixed-based 
biochar on soil water infiltration 


Cumulative soil water infiltration refers to the total amount of water infiltrating into the soil 
through a surface unit within a specific time after the beginning of infiltration, and it is the 
integral of the infiltration rate as a function of time. Before the soil water infiltration process is 
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stable, cumulative soil water infiltration is often used to characterize soil water infiltration 
capacity. Figure 4 shows that cumulative soil water infiltration of each treatment gradually 
increased with time, but the effects of treatments on cumulative soil water infiltration were 
different. Compared with CK, all other treatments reduced cumulative soil water infiltration, 
except for B4%s2). For mixed-based biochar with S1 particle size, with the increase of 
application amount, the inhibitory effect on soil water infiltration first strengthened and then 
weakened under application patterns A and C. Mixed-based biochar with S2 particle size had the 
same effect on soil water infiltration inhibition, which was first unchanged and then increased. 
Cumulative soil water infiltration increased in the order of 2% s1)<4% s2)< 1% s1)<4%s1)< 
1% (s2)<2% s2)<CK. For mixed-based biochar with S1 particle size, inhibition of soil water 
infiltration under application pattern B first weakened and then remained constant gradually with 
the increase of application amount; however, for mixed-based biochar with S2 particle size, with 
the increase of application amount, inhibition of soil water infiltration first remained constant and 
then weakened sharply. Therefore, the treatment of 1% s1) showed the best inhibitory effect, and 
the treatment of 4% ,s2) did not inhibit water migration but promoted soil water infiltration. The 
remaining treatments had little effects soil water infiltration, and the effects increased in the order 
of 1% s1)< 1% (s2)<2% (82) and 4% (51)<2%(s1)<CK<4% g2). 


— CK — 1%s1) —-s 1%s2) —- 2%Q1) 2% s2) n 4%s1) 4% 52) 


(a) Application pattern A 20 (b) Application pattern B 


(c) Application pattern C 
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infiltration (cm) 
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Fig. 4 Effects of application patterns, application amounts, and particle sizes of mixed-based biochar on 
cumulative soil water infiltration. (a), application pattern A; (b), application pattern B; (c), application pattern C. 


With the same application amount, the effects of different application patterns on cumulative 
soil water infiltration were different. When application amount was 1%, the order of the effects of 
application patterns on cumulative soil water infiltration was A<B<C, and S1 particle size 
reduced cumulative soil water infiltration. In short, for any application amounts, cumulative soil 
water infiltration under application pattern A was lower than those under application patterns B 
and C, and cumulative soil water infiltration under S1 particle size was lower than that under S2 
particle size. 

Considering the influences of application patterns, application amount, and particle sizes of 
mixed-based biochar on cumulative soil water infiltration, an infiltration time of 10 min was 
taken as an example. Among all the treatments, only the treatment of B4%ç(s2) showed a lower 
infiltration (13.5 cm), which was greater than the CK (10.4 cm) (P<0.05), and the increase in 
cumulative soil water infiltration was 30% for CK, which promoted cumulative water infiltration 
of aeolian sand soil. Cumulative soil water infiltration under different treatments varied as 
follows: A2% (s1)<A4%(s1)< Al % (51)< A 4% s2)<C2%(51)<B 1 %(s1)< Al % (52)<C4.%(s2)<B 1 %(52)< 
C4%(§1)< C1 %(51)< A2% (52) < C1 %(52)< B2 %(g2)<B4-%(51)<B2%(51)<C2%(s2)<CK. Among them, 
treatments of A2%s1), A4%s1), A1%s1), A4% 52), C2% 51), and B1%,s1) showed the best effects, 
and their reduction effects were all higher than 30%, which met the requirements of loess soil 
hydraulic parameters suitable for plant growth. 
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3.3 Effects of application patterns, application amounts, and particle sizes of mixed-based 
biochar on soil water infiltration parameters 


Several theory-based approximate and empirical models that simplify the concepts involved in 
infiltration have been developed for field applications to further study the effects of application 
patterns, application amounts, and particle sizes of mixed-based biochar on soil water infiltration 
(Table 5). 


Table 5 Models used to determine the effects of mixed-based biochar on soil water infiltration (Wang et al., 
2017) 


Model name Equation Model parameter 


I is the cumulative soil water infiltration (cm); K, is 

the saturated hydraulic conductivity (cm/s); 6; and 0, 

Green and Ampt (1911) 7=K,01+(@, -6)S; Ir] are the presumed initial and saturated water contents 

Theory-based (g/g), respectively; and Sy is the average potential 
model suction of the wetting front (cm). 


S is the sorptivity (cm/s®°); A is the stable 


Philip (1957 = 519° 
ilip-(1957) FEST RAL permeability (cm/s); and ¢ is the infiltration time (s). 
: K (K>0) and n (O<n<1) are the dimensionless 
Kostiakov (1932) T=Kt" empirical constants. 
a and b are the presumed initial and final soil water 
i ct infiltration rates (cm/s), respectively; c is an 
Empirical Horton (1941) I=at+1/c(b-a)\(l-e ©) empirical constant; and e is the natural logarithm 
model (approximately 2.71828). 
USDA-NRCS (1974) I =a"t?" +0.6985 a" and b" are the dimensionless empirical constants. 


K' (K'>0) and b' (0<b'<1) are the dimensionless 
1=K't’ +A't empirical constants, and A' (A'>0) is the final soil 
water infiltration rate at the steady state condition. 


Kostiakov-Lewis 
(Mezencev, 1948) 


Note: The USDA-NRCS model is the United States Department of Agriculture-Natural Resources Conservation Service model. 


According to Liu et al. (2010), an implicit functional relationship between cumulative soil 
water infiltration and infiltration time exists in the Green and Ampt model and it should be 
reflected in the infiltration equation. However, this process may introduce new parameters into 
the adjustment equation and cause errors, therefore, the Philip model was the only theory-based 
model selected for this study. 

The characteristics of soil water infiltration under different application patterns, application 
amounts, and particle sizes of mixed-based biochar were compared and analyzed. Table 6 lists the 
estimated parameter values for the five infiltration models. It was concluded that there were 
different infiltration modes in aeolian sand soil for the five infiltration models with the same 
experimental treatments. The values of K (saturated hydraulic conductivity (cm/s)), b (presumed 
final infiltration rate (cm/s)), a" (dimensionless empirical constant), and K' (dimensionless 
empirical constant) were very similar. These empirical parameters had similar soil physical 
meanings with that of the sorptivity for Philip's theoretical basic model (adsorption). A 
comparison of the parameter values for each model showed that K and K' were more significant 
than a", which was similar to the results reported by Wang et al. (2017). 

Table 6 shows that for all treatments, 95% of the R? values exceeded 0.900, indicating perfect 
performance of the five models. The Horton model was the best, with R? exceeding 0.990 (except 
for the treatment of B2%,s2)), which was similar to the results reported by Duan et al. (2011). At 
the same time, we found that the fitting degrees of the five models under S1 particle size were 
higher than those under S2 particle size, and those under application pattern C were significantly 
better than those under application patterns A and B. 

Taking application pattern C as an example, the parameters generated by each infiltration 
model for the different treatments were similar to the K values of the Kostiakov model. According 
to the fitting results, the K values of biochar treatments were less than that of CK, indicating that 
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application amounts and particle sizes of mixed-based biochar were beneficial for inhibiting 
initial soil water infiltration. However, as application amount of mixed-based biochar increased, 
the K values first decreased and then increased, indicating that excessive application of 
mixed-based biochar showed a negative effect on soil water infiltration, reducing the water 
infiltration capacity of aeolian sandy soil. Moreover, the K values for different particle size 
treatments did not show a consistent pattern. Additionally, under application pattern C, by fitting 
the measured data of each treatment with the Philip infiltration model, it was found that the 
absorption rate and the stable soil water infiltration rate of different application amounts and 
particle sizes were significantly different from those under CK. This showed that application 
amount and particle size of mixed-based biochar all affected the process of soil water infiltration. 
The absorption rates under all biochar treatments were lower than that under CK, indicating that 
biochar treatments significantly reduced soil capillary suction and their abilities to absorb or 
release liquid were also weaker than that of CK. With the increases of application amount of 
mixed-based biochar, the influence laws for different particle sizes were not consistent. Further, 
application amount and particle size significantly increased soil water infiltration rate and 
improved water infiltration performance of aeolian sandy soil. Specially, the treatments of 
C1% s2) and C2%,s2) were the most effective among all treatments. 


Table 6 Parameters in the five infiltration models 


Kostiakov-Lewis 


Philip Kostiakov Horton USDA-NRCS 
(Mezencey) 


Treat- 
ment 


A R? K n R? a b c R? a" b" R? K' A' b' R? 


CK 2.972 0.163 0.919 1.837 0.862 0.868 0.045 2.935 0.280 0.993 2.672 0.562 0.936 3.629 0.010 0.453 0.940 
Al%s1) 1.481 0.145 0.976 1.145 0.730 0.984 0.110 1.075 0.129 0.999 1.052 0.711 0.975 1.455 0.005 0.625 0.986 
A1%s2) 1.684 0.149 0.969 1.164 0.760 0.973 0.081 1.112 0.103 0.997 1.236 0.693 0.971 1.507 0.003 0.650 0.980 
A2% s1) 1.621 0.042 0.957 1.072 0.700 0.971 0.082 0.968 0.125 0.995 1.096 0.639 0.956 1.449 0.010 0.574 0.969 
A2% s2) 2.148 0.010 0.928 1.683 0.600 0.934 0.078 1.465 0.151 0.994 2.187 0.476 0.933 2.527 0.001 0.453 0.942 

0.090 0.999 1.301 0.621 0.964 1.682 0.003 0.564 0.974 


A4% s2) 1.532 0.010 0.866 1.607 0.510 0.910 0.046 1.123 0.135 1.000 2.081 0.399 0.935 3.101 0.002 0.310 0.932 


A4% s1) 1.836 0.043 0.964 1.084 0.720 0.968 0.053 0.966 


Bl1%(sı) 1.092 0.343 0.989 1.089 0.810 0.994 0.110 1.039 0.093 0.992 0.895 0.839 0.985 1.091 0.150 0.706 0.992 
B1%s2) 1.400 0.292 0.986 1.131 0.810 0.992 0.150 1.056 0.078 0.999 1.115 0.776 0.987 1.446 0.020 0.701 0.992 
B2%s1) 1.429 0.430 0.967 1.168 0.870 0.992 0.287 1.121 0.061 0.996 1.280 0.805 0.972 1.254 0.051 0.814 0.992 
B2%s2) 2.267 0.010 0.915 1.357 0.700 0.939 0.010 1.459 0.111 0.981 2.126 0.519 0.908 2.200 0.001 0.528 0.913 
B4% si) 1.925 0.267 0.974 1.290 0.810 0.988 0.032 1.167 0.054 0.998 1.532 0.721 0.979 1.731 0.003 0.710 0.978 
B4% 52) 2.392 0.010 0.483 3.185 0.450 0.736 0.029 2.775 0.185 0.996 5.660 0.266 0.829 8.513 0.010 0.158 0.960 
C1%s1) 0.735 0.551 0.996 1.066 0.879 0.997 0.051 1.062 0.067 1.000 0.790 0.955 0.987 1.173 0.001 0.828 0.998 
C1%s2) 0.448 0.728 1.000 1.145 0.875 1.000 0.820 1.410 1.044 1.000 0.703 1.061 0.991 0.445 0.768 0.386 1.000 
C2%s1) 1.070 0.327 0.994 1.075 0.803 0.992 0.310 1.088 0.167 1.000 0.857 0.841 0.988 1.271 0.102 0.652 0.997 
C2%s2) 0.909 0.609 0.984 0.989 0.958 0.992 0.066 1.144 0.042 0.996 1.044 0.900 0.984 1.298 0.021 0.838 0.989 
C4%s1) 1.100 0.396 0.992 1.127 0.827 0.996 0.089 1.069 0.074 1.000 0.959 0.847 0.989 1.323 0.012 0.752 0.996 
C4%s2) 1.530 0.228 0.984 1.237 0.760 0.990 0.092 1.115 0.093 1.000 1.169 0.739 0.984 1.533 0.001 0.668 0.991 


3.4 Effects of application patterns, application amounts, and particle sizes of mixed-based 
biochar on soil water evaporation 


Surface water evaporation is a significant cause of soil water loss. Effectively inhibiting 
evaporation of soil water is of great significance in improving soil water use efficiency. Figure 5 
shows the variation trend of cumulative soil water evaporation with different evaporation times. 
Cumulative soil water evaporation under the treatment of B4%s1) was similar to that under 
CK. The remaining treatments all reduced cumulative soil water evaporation. Under application 
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patterns A and B, a significant difference appeared after 25 d of continuous soil water 
evaporation. Specifically, under S1 particle size, cumulative soil water evaporation first decreased 
and then increased with increasing application amount, while that under S2 particle size was 
constant at first and then increased. Cumulative soil water evaporation increased as 2% s2)< 1% s2), 
4% (52), 2% s1), and 1% s1)<4%s1) and CK. Under application pattern C, there was a significant 
difference after 20 d of continuous soil water evaporation, and cumulative soil water evaporation 
was inversely proportional to application amount and particle size of mixed-based biochar. 
Cumulative soil water evaporation increased as 4% s2)<4% s1)<2.%(s2)<2%(s1)< 1% s2)< 1 MH s1)<CK. 
Compared with CK, the treatment of 4% s2) significantly reduced cumulative soil water 
evaporation by 33%. For the same application amount, the difference of cumulative soil water 
evaporation between application patterns A and B was slight. Cumulative soil water evaporation 
rates under application patterns A and B were significantly larger than that under application 
pattern C. Under application patterns A, B, and C, cumulative soil water evaporation rates of 
mixed-based biochar with S1 particle size were smaller than those of mixed-based biochar with 
S2 particle size. Mixed-based biochar reduced ineffective soil water evaporation in the field and 
improved the water holding capacity of aeolian sandy soil. The combination of application 
pattern, application amount, and particle size was still valuable (P<0.05), which had almost no 
effect on cumulative soil water evaporation; specifically, cumulative soil water evaporation 
increased as C4%(82)<C4-% 1g 1)<C2%(82)< C2%(g1)<A2%(s2)<C 1 %(§2)<B2%(82)<C 1 Gg 1)< A4-%s2) 
<A1% (s2)< A1%(s1)<B1 %(§2)< A2%(81)<B4-%(82)<B 1 Pens A4-%(51)<B2%s1)<CK, C4%s1), C4% 52), 
C2% 81) and C2% 82). 
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Fig. 5 Effects of different application patterns, application amounts, and particle sizes of mixed-based biochar 
on cumulative soil water evaporation 


4 Discussion 


4.1 Main and interactive effects of application patterns, application amounts, and particle 
sizes of mixed-based biochar on soil water infiltration and evaporation 


In this study, we designed a mixed-level orthogonal experiment to explore the effects of three 
influencing factors of mixed-based biochar (application pattern, application amount, and particle 
size) on water infiltration and evaporation of aeolian sandy soil. According to a visual analysis of 
the mixed-level orthogonal experiment, we analyzed the order and interactions for the effects of 
the three influencing factors. It was found that the order for the effects of influencing factors on 
soil water infiltration decreased as application pattern>application amount>particle size. This 
process involved the following interactions: application amountxparticle size>application 
patternxparticle size>application patternxapplication amount. However, there was no interaction 
of the effects of influencing factors on soil water evaporation, and the single effect decreased as 
application pattern>particle size>application amount. Therefore, application pattern of 
mixed-based biochar was the main factor affecting soil water infiltration and evaporation. The 


ZOU Yiping et al.: Effects of mixed-based biochar on water infiltration and evaporation... 


reason was that the soil-carbon interface directly affected the content, hydraulic gradient, and 
conductivity of soil water. Table 6 shows that for the whole infiltration process, the difference 
between absorption rate and stable soil water infiltration rate for application patterns were more 
significant than those for application amounts and particle sizes, which also verified the 
above-mentioned conclusions. Moreover, the changes of steady soil water infiltration rates were 
more significant than those of absorption rates, implying that the effect of mixed-based biochar on 
stable soil water infiltration was greater in the later stage than in the initial stage. Additionally, 
application amount and particle size directly affected soil porosity, so the interaction was 
significant. In this study, application amount changed the soil porosity more significantly, which 
depended on the application amount and particle size of mixed-based biochar. Whether there was 
a synergistic or antagonistic effect between the two factors also depended on the levels of them. 
Figure 5 shows that in the evaporation process, the application of mixed-based biochar had a 
slightly more significant impact in the first stage than in the second stage. Different application 
patterns, especially application pattern C, significantly affected soil water evaporation in the two 
stages. The effect of particle size on soil water evaporation was greater than that of application 
amount, because particle size had a more direct and significant impact on capillary action. 


4.2 Effects of application patterns, application amounts, and particle sizes of mixed-based 
biochar on soil water infiltration 


The soil water infiltration process was mainly affected by the water supply intensity (external 
factors) and soil water infiltration capacity (soil texture, structure, bulk density, and initial 
moisture content). This water supply intensity and the initial moisture content were consistent in 
the study. The soil water infiltration capacity mainly depended on the porosity of the soil. After 
mixed-based biochar with a rich pore structure and high specific surface area was applied to the 
soil, it first changed soil bulk density, increased soil porosity, influenced soil aggregates and their 
stabilities, and ultimately affected the hydraulic properties of the soil (Xiao et al., 2015). 

When a lower amount of 1% mixed-based biochar was added, it had a more significant impact 
on soil texture, bulk density, and porosity. Moreover, mixed-based biochar with finer particle 
sizes filled some of the soil pores. Compared with pure soil, mixed-based biochar with rich pore 
structure and hydrophilic characteristics of inhibited the downward migration of the water 
(Nguyen et al., 2010). Therefore, the inhibitory effect of application pattern A (0—20 cm) on water 
infiltration was more evident than those of application pattern B (0-10 cm) and application 
pattern C (10-20 cm). When application amount was as low as 1%, the inhibitory effect was not 
obvious due to the small amount of mixed-based biochar. Both the carbon-soil mixed layer and 
the pure soil layer had essential impacts on soil water infiltration. For application pattern B, the 
moisture absorption capacity of the surface layer was weaker than that of the underlying layer. 
Because of the inhibiting effect of the surface layer on soil water infiltration, the speed of soil 
water infiltration to the carbon-soil interface was slow, which triggered advantageous moisture 
absorption by the pure soil in the underlying layer. Under application pattern C, when the suction 
reached the wetting front and the carbon-soil interface, the underlying layer suction was less than 
the surface layer suction and the matric suction at the wetting front could not meet the 
requirement for rapid soil water infiltration. However, the matric suction at the wetting front 
decreased continuously with the supply of water. When the matric suction reached the mixed soil 
layer, the wetting front entered the mixed soil layer (Wang et al., 2010). However, the moisture 
content of the surface soil was higher under application pattern C, which increased the infiltration 
speed of the wetting front to a certain extent through the carbon-soil interface and made the 
overall infiltration speed faster than that under application pattern B. When medium-high 
application amounts (2% and 4%) of mixed-based biochar were applied, the carbon-soil mixed 
layer significantly inhibited water infiltration, and the influence of the pure soil on water 
infiltration was relatively weakened. Therefore, when the soil layer thickness remained the same, 
the effects of application patterns B and C gradually approached each other. 

Under application patterns A and C, the distribution characteristics of the soil pores were 
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changed as the initial addition of mixed-based biochar with S1 particle size. The total volume of 
macropores was reduced, while the total volume of small and medium pores was increased. 
Moreover, the water channel was tortuous and complex, which significantly inhibited soil water 
infiltration (Wang et al., 2015). When application amount of mixed-based biochar continued to 
increase up to 4%, some biochar stuck together or adsorbed on the soil and formed tiny pellets 
because of the sharp increase in the amount of mixed-based biochar. Compared with 2% 
application amount, the volume of macropores was relatively increased under 4% application 
amount, and the inhibition of soil water infiltration was no longer significant. For mixed-based 
biochar with S2 particle size, the reductions in soil macropore volumes were similar under 1% 
and 2% application amounts, which may be related to the finer particle size and the slight change 
in application amount. Further, when application amount of 4% mixed-based biochar was added, 
some of the soil pores were filled by mixed-based biochar with finer particle size. Even if 
mixed-based biochar with S2 particle size was bonded, it could not form larger aggregates. In 
addition, the rich pore structure and hydrophilic characteristics of mixed-based biochar itself 
could more effectively inhibit the infiltration of soil water (Wen and Zheng, 2012). However, the 
overall inhibitory effect of application pattern A on soil water infiltration was more potent than 
that of application pattern C, mainly because application pattern C contained pure soil in the 0—10 
cm surface layer. The volume of large pores in the soil was larger than that in the carbon-soil 
mixture, so the pure soil in the 0-10 cm surface layer showed poor suppression of soil water 
infiltration. In addition, for cultivated soil, there was a sizeable gap in the matrix potentials in the 
soil-carbon interface, which promoted the absorption of the carbon-soil mixture under application 
pattern C and ultimately inhibited soil water infiltration in the 0-20 cm layer (Li, 2019). With the 
increase in application amount of mixed-based biochar with S1 particle size under application 
pattern B, the inhibition of soil water infiltration first sharply weakened and then gradually 
remained unchanged. In contrast, with the increase in application amount of mixed-based biochar 
with S2 particle size under application pattern B, the soil water infiltration first unchanged and 
then weakened sharply. Soil water infiltration was affected by the soil water infiltration ability 
and water supply intensity. In the initial stage of soil water infiltration, when the water supply of 
the Mahalanobis flask was unstable, the water supply intensity was extremely high. The matric 
suction at the wetting front could not meet the requirements for rapid soil water infiltration, 
causing water outflow from the surface layer of the soil column. When the application of 
mixed-based biochar was increased sharply, soil water infiltration was more difficult, resulting in 
water overflow. This led to the greatest suppression at the low application amount (1%) for 
mixed-based biochar with S1 particle size. For mixed-based biochar with S2 particle size, the 4% 
application amount caused an increase in cumulative soil water infiltration due to water overflow, 
while cumulative soil water infiltration rates were similar for other treatments. 


4.3 Effects of application patterns, application amounts, and particle sizes of mixed-based 
biochar on soil water evaporation 


Soil water evaporation is the main loss caused by hydrological cycles in arid and semi-arid 
regions. Soil properties such as water content, texture, structure, and color all affect soil water 
evaporation. 

For any application amount of mixed-based biochar, we found that the evaporation process 
could be divided into two stages. Compared with application patterns A and B (25 d), soil water 
evaporation reached stability earlier (20 d) under application pattern C, and evaporation did not 
increase sharply. Under application pattern C, soil water evaporation rate in the first stage (the 
beginning of soil water evaporation) was lower than those under application patterns A and B. In 
the first stage, water was first quickly lost from the larger pores, so soil water evaporation 
completed earlier under application pattern C than under application patterns A and B because of 
the large pores in the surface layer. In addition, soil water evaporation entered the second stage 
earlier under application pattern C and the soil water evaporation loss in this stage occurred 
through steam diffusion. Since the soil-carbon interface blocked the capillary pores, the pore 
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connectivity of the whole soil decreased, and it was difficult to increase the diffusion rate of water 
vapor to reduce soil water evaporation. Therefore, covering the surface with coarse aeolian sandy 
soil could help to reduce soil water evaporation, which was consistent with the research of Wang 
et al. (2018). Moreover, the complex internal pore structures developed in the surface layer of 
mixed-based biochar with application patterns A and B increased the number of capillary pores in 
the soil, which reduced the number of large soil pores, relative to that of the pure soil in the 
surface layer of mixed-based biochar with application pattern C. Xia et al. (2000) found that soils 
with fine texture and well-developed capillary pores can promote the movement of water from the 
underlying layer to the surface layer. Evaporation pulls water during the evaporation process (Xia 
et al., 2000). Hence, the effects of application patterns A and B on soil water evaporation were 
more lasting and more robust than that of application pattern C. For mixed-based biochar with S2 
particle size, cumulative soil water evaporation under application patterns A and B in the second 
stage decreased first and then increased with the increase of application amount, while for 
mixed-based biochar with S1 particle size, cumulative soil water evaporation appeared to be 
unchanged first and then increased. In the second stage, evaporative loss occurred through 
capillary pores. With a high application amount of mixed-based biochar, the capillary pores of the 
soil increased significantly in the second stage. Finally, cumulative soil water evaporation under 
4% application amount of mixed-based biochar was gradually close to that under CK. As more 
mixed-based biochar was applied, the finer the particle size was, the larger the specific surface 
area, the stronger the water adsorption, and the more the stored water. When the interface blocked 
the capillary pores and the pore connectivity of the entire soil decreased, the water stored in the 
underlying layer could not migrate to the surface layer and evaporate upward, which finally 
reduced soil water evaporation (Li et al., 2016). 


5 Conclusions 


The soil used in this study was aeolian sandy soil from a mining area. With heavy precipitation, 
soil erosion was prone to occur, which increased the rate for the loss of soil nutrients. This has 
affected the growth of crops and exacerbated the fragile ecological environment. Appropriate 
application of mixed-based biochar can improve soil water infiltration in the mining areas and 
meet the requirements of loess soil hydraulic parameters suitable for plant growth. We separately 
used five infiltration models to fit cumulative infiltration and time, and found that although the 
Horton model fitted the effects of application pattern C and S1 particle size best, it was suitable 
for almost all treatments. Additionally, we found that applying mixed-based biochar to the 
underlying layer (10-20 cm) can effectively prevent soil water evaporation. After applying 2% of 
mixed-based biochar with S1 particle size to the underlying layer, cumulative soil water 
infiltration and evaporation were significantly improved while minimizing water evaporation. 
This enhanced the water retention capacity of aeolian sandy soil and enabled the efficient use of 
limited water resources. Moreover, we found that the application pattern was the main factor 
affecting soil water infiltration and evaporation. There were interactions among factors 
(application pattern, application amount, and particle size) in the infiltration process, and the 
interaction between application amount and particle size was the strongest. In contrast, there was 
no interaction among these three factors for the evaporation process. However, this conclusion is 
based on a short-term laboratory column experiment. Therefore, it is necessary to further study 
the effects of long-term application of mixed-based biochar on soil physical properties and plant 
growth in mining areas. 
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